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Damped Population Oscillation in a Spontaneously Decaying 
Two-Level Atom Coupled to a Monochromatic Field 
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We investigate the time evolution of atomic population in a two-level atom driven by a monochro- 
matic radiation field, taking spontaneous emission into account. The Rabi oscillation exhibits am- 
plitude damping in time caused by spontaneous emission. We show that the semiclassical master 
equation leads in general to an overestimation of the damping rate and that a correct quantitative 
description of the damped Rabi oscillation can thus be obtained only with a full quantum mechanical 
theory. 

PACS numbers: pacs 



The atom-field interaction is one of the most funda- 
mental problems of quantum optics El • Despite 
the success of the semiclassical approach which treats 
the atom quantum mechanically but the field classically, 
there exist many optical phenomena, an accurate descrip- 
tion of which requires a full quantum-mechanical treat- 
ment of both the atom and the field. One such example 
is the collapse and revival 0, 0, El> llfl of the Rabi 
oscillation in the Jaynes-Cummings model 

The influence of dissipation (spontaneous decay, cav- 
ity damping ) on the collapse and revival of the Rabi 
oscillation has been studied in the past 0, 0] . When 
one considers a transition in the microwave region where 
the first observation of the collapse and revival was made 
0|, spontaneous decay can certainly be neglected. The 
past investigations on the influence of the cavity damp- 
ing on the collapse and revival have shown that the col- 
lapse is not affected much by the cavity damping, but 
the revivals, especially at large times, suffer attenuation 
in the oscillation amplitude [12L . In the optical do- 
main, spontaneous emission may be expected to play a 
much stronger role. It has been shown, however, that re- 
vivals in the optical region are still much more sensitive 
to cavity decay than to spontaneous emission (lif. 

In this letter we report on our study of the influ- 
ence of spontaneous emission on the collapse and re- 
vival of the Rabi oscillation. The reason for our study 
of spontaneous emission, despite its relative unimpor- 
tance compared with cavity decay, is that there exists 
a subtle quantum effect when the collapse and revival 
phenomenon is combined with spontaneous decay. The 
quantum effect has to do with the damping rate of the 
amplitude of the Rabi oscillation. As is well known, the 
amplitude of the Rabi oscillation decreases progressively 
at each successive revival, even in the absence of sponta- 
neous decay and any other dissipation, because a smaller 
and smaller number of dephased oscillations rephase at 
each successive revival time. Spontaneous decay works 
to further decrease the oscillation amplitude over the al- 
ready existing decrease in its absence. We are in this 
work mainly concerned with the damping rate of the 



oscillation amplitude due to spontaneous emission. We 
show that the semiclassical theory tends to overestimate 
this damping rate. A correct quantitative description of 
the damped Rabi oscillation is thus obtained in general 
only with a full quantum mechanical treatment. The 
overestimation stems from the fact that the semiclassi- 
cal picture of the transitions involved is oversimplified, 
as will be seen below. 

Let us first look at the semiclassical treatment of a two- 
level atom driven by a resonant, monochromatic radia- 
tion field. The driving field induces the Rabi oscillation 
between the excited state |e) and the ground state \g) of 
the atom. We also take into consideration spontaneous 
decay of the excited state |e) into the ground state \g) 
and other external states. The transitions that the atom 
undergoes are schematically shown in Fig. ^ The master 
equation for the density operator p of the atom is given 
in the interaction picture under the dipole and rotating 
wave approximations by 

dp 

dt 

r 

+ o(26cr_/?(T+ - a+a-p- pcr+CT-), (1) 



— = -i[na+ + fTcr_,p] 



where fl is the Rabi frequency, <r + and er_ are the 
atomic raising and lowering operators given respectively 
by <r + = | e) (<7 1 and cr_ = |<?)(e|, T represents the total 
spontaneous decay rate of the excited state, and the pa- 
rameter b (0 < b < 1) is the branching ratio into the 
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FIG. 1: Stimulated and spontaneous transitions of a two-level 
atom in the semiclassical picture, fl is the Rabi frequency and 
bT and (1 — b)T represent the spontaneous emission rate into 
the ground state and other external states, respectively. 
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ground state, i.e., the decay rate into the ground state is 
bT and that into all other external state is (1 — b)T. 

From Eq. we immediately obtain the optical Bloch 
equations for the density matrix elements p gg = {g\p\g), 
p ee = (e\p\e), and p ge = (g\p\e), 



Pgg 

Pee 
Pge 



-i (fi* peg - n Pge) + bTpee, (2a) 
-i (fi p g e - ST Peg) ~ Tpee, (2b) 



-in*{pee - Pgg) ~ -Pge 



(2c) 



Eqs. J2J) yield a third-order differential equation for p e 
which can be solved analytically to give 



p ee (t) = e^V 1 '- 
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(4a) 
(4b) 



a = 216 &r|nr + 3^/5184 6 2 r 2 |^r + 3 (16 |n| 2 -r 2 )- 

(4c) 

In obtaining Eq. |JSJ, we have assumed that the atom 
is initially in its ground state; p gg (0) = 1 and p ee (0) = 

Pge(0) = 0. 

From now on we limit ourselves to the strong cou- 
pling regime in which fi is sufficiently greater than T 
and, as a result, both fa and fa are real and positive. 
Eq. © indicates that, in the strong coupling regime, the 
excited state population exhibits a damped oscillatory 
behavior with a frequency of ^fa and a damping rate 
of (r + /i)/2. Our calculation indicates that in general 
fa increases roughly linearly with b while fa does not 
depend much on b, which means that the damping rate 
of the population oscillation increases with b while the 
oscillation frequency is roughly independent of b. 

As a special case we set b = in Eq. ® and obtain 
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(5) 



where £ = y4|fi| 2 — In the other extreme case in 
which b = 1 , we obtain 

4|fi| 2 r „ , , sr . 3ptl , . 
Pee{t) = 8io p + r 2 [ " (cos xt + 4A } ] ' (6) 



becomes 



where, A = \/4|0| 2 - In the limit |0| > T, Eq. © 



p ee (i) £ -(l-cos(2|0|t)e-i rt ). 



(7) 



The solutions, Eqs. (J5J and (JHJ, have been given previ- 
ously |3, although, to our knowledge, the general so- 
lution, Eq. |J3J), has not. One sees from Eqs.JSJhQ that 
the semiclassical treatment predicts that the damping 
rate of the population oscillation increases with respect 
to the branching ratio b from ^ at b — to |r at b = 1. 

Why does the population oscillation damp out faster 
when the branching ratio b is greater, i.e., when the spon- 
taneous decay to the ground state has a higher weight? 
The clue to this question can be provided by examin- 
ing the role played by the last term bTp ee of Eq. 
Roughly speaking, this term represents a kind of pump- 
ing of population into the ground state. This "inter- 
nal pumping" into the ground state, being proportional 
to pee, is strongest when the excited state population 
is maximum, i.e., when the ground state population is 
minimum. The transition from the excited state to the 
ground state occurs because the excited state population 
is greater, but this "internal pumping" due to the term 
bT pee works in the direction to spoil the very reason for 
the transition. As a result, the greater this term bTp ee 
is, the faster the population oscillation dies out. 

We now present a full quantum-mechanical description 
treating the driving field as well as the atom quantum me- 
chanically. The master equation for the density operator 
p of the atom-field system reads 

dp . r . * t l 

— = -i[ga + a + g cr_a',pj 

r 

+ 2 (26cr_/?cr + - (j+o-p - pv+cr-), (8) 

where and a are, respectively, the creation and anni- 
hilation operators of a photon of the driving field mode, 
and g is the parameter that measures the atom-field cou- 
pling (not to be confused with the ground state). 
The optical Bloch equations now take the form 



Pgn+l,gn-\-l 
Pen, en 
Pgn-\-l,en 



-iVn + 1 (g*p en ,gn+l - 9 Pgn+l,en) 
+ bTp en+ ien+l, (9a) 
—iVn + 1 (g p gn +l,en - g* Pen,gn+l) 
T Pen, en i 

(9b) 



-Wn + lg*(p en ,en - Pgn+\,gn+l) 

r 

7\ Pgn-\-l,ent 

(9c) 

where p gn +i,gn+i = (.9, n + 1| p \g, n + 1), p en ,en = 
(e,n\p\e,n), p gn+ \, en = (g,n + l|p|e,n), and \e,n) rep- 
resents the state of the atom-field system in which the 
atom is in the excited state and the field has n photons, 
and similarly for \g, n + 1). 

The transitions that are represented by Eqs. © are 
depicted in Fig. [21 We emphasize that the transition be- 
tween \g,n+ 1) and \e,n) is not closed even for b = 1, 
because spontaneous emission of a photon from \e, n + 1) 
increases the population of the state \g,n+ 1), and simi- 
larly spontaneous emission of a photon into \g, n) as well 



FIG. 2: Stimulated and spontaneous transitions of the atom- 
field system in the quantum picture. Only the transitions 
that directly influence the populations of the states |e, n) and 
\g,n + 1) are shown. |gr| is the parameter that measures the 
atom-field coupling, and &r and (1 — b)T represent the sponta- 
neous emission rate into the ground state and other external 
states, respectively. 



FIG. 3: Variation of the excited state population with 
time, (a) p%% act (t) (green-dotted), pZ ml {t) (red-solid), and 
pi" ant {t) (blue-dashed) as a function of time for the case 
\g\ = 10F, (n) = H 2 = 25 and b = 1. ptT (t) and p q e u e ant {t) 
are indistinguishable, (b) an expanded view of the revival 
region. 



as into external states decreases the population of the 
state |e,n). This "openness" of the quantum transitions 
is responsible for the quantum effect we discuss in the 
following. 

Eqs. © cannot be solved analytically in general, be- 
cause, as mentioned above, the transition between \g, n + 
1) and |e,n) is not completely closed. Only in the limit 
b = 0, Eqs. © are reduced to a closed set of equations 
relating p s „+i, 9 „+i, /o e „, e „, and p gn +i,en, and an ana- 
lytic solution for p ee {t) = S^Lo Pen,en{t) can be obtained 
which reads 



Pee(t) = 



4|<?| 2 (n + l)e-l tt l 2 a 2 (™ +1 ) 



_E. . o £n , 



71=0 



(n + l)!|£„|2 ™ 2 

(10) 

where £ n = ^ A\g\ 2 (n + 1) - In obtaining Eq. (JTJU, 
we have assumed that initially the atom is in the ground 
state and the driving field in the coherent state |a), i.e., 

_ j a |2 2(71 + 1) 

Pgn+l, 97l+l(0) = " j Pen, en (0) = Pgn+l,en{0) = 

0. Comparing Eq. (|10|) with its semiclassical counter- 
part, Eq. (J5J, we first note that Eq. UlUf) consists of the 
summation of an infinitely large number of terms arising 
from the fact that the Rabi frequency for the transition 
between \g,n+ 1) and |e, n) is dependent on n. As is 
well known, this leads to the collapse and revival of the 
Rabi oscillation |E IE IE IE IE E3 • Our main concern 
in this work is however not the collapse and revival but 
the damping rate of the population oscillation. Eq. l(TT)|) 
indicates that this damping rate for the case b = is 
2, in agreement with the semiclassical damping rate of 
Eq. (J5J. The question then arises: does the damping rate 
calculated according to the quantum theory agree with 
the semiclassical damping rate also for nonzero values of 
&?, i.e., does the quantum damping rate increase with the 
branching ratio b as predicted by the semiclassical the- 
ory? This is the main issue we wish to address ourselves 
in this work. 

We show in the following that the answer to the above 
question is no and that the quantum damping rate of 



the population oscillation does not vary much with the 
branching ratio b. In order to see this, let us first re- 
call that the reason for the increase of the damping rate 
with respect to the branching ratio in the semiclassical 
theory is that the "internal pumping" represented by the 
last term bTp ee of Eq. (0i) works to weaken the oscilla- 
tion. Similarly, we see from Eq. J^t) that the last term 
kTpen+i, e7i+i P ia y s the r °l e of pumping population into 
the state \g,n + 1). The strength of this "pumping" is 
determined by the population of the state \e, n+ 1), while 
the state \g, n + 1) into which population is pumped un- 
dergoes the Rabi oscillation with the state |e, n) not with 
the state \e,n + 1). Since there is no definite phase re- 
lationship between the population of the state |e, n + 1) 
and that of the state |e,n)(or \g,n + 1)), this pumping 
strengthens the oscillation as much as it weakens the os- 
cillation. After a sufficiently long time, therefore, the 
effect of the pumping term can largely be neglected. We 
thus expect that the damping rate of the population os- 
cillation at any nonzero value of b is roughly the same as 
that at b = 0, i.e., it is £ regardless of b. This expecta- 
tion based on the quantum-mechanical consideration of 
transitions contrasts sharply with the semiclassical pre- 
diction. 

We now present numerical data to support the above 
expectation. The "exact" solution for p ee (t), which we 
denote by p e e ^, act {t), was obtained by numerically inte- 
grating the optical Bloch equations, Eqs. ©. The system 
parameters we chose for our computation are \g\ = 10r, 



= 25, and 6=1. The ratio 



10 can be 



achieved in an optical cavity [1E0. The value b = 1 
is chosen because at this value the semiclassical damp- 
ing rate takes on a maximum value of |r and differs 
most from the expected quantum damping rate of £. In 
Fig. G2a) we plot p^ act as a fuction of time. We try to 
fit this exact curve with an approximate formula using a 
guessed value T guess of the damping rate. We take the 
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FIG. 4: The differences S semi 
dotted) and 8 quant (t) = pZ act 
the case \g\ = 10r, (n) = \a\ 2 



(t) = p%T c \t) -p««™(*)(red- 
(t) - pie ant (t) (blue-solid) for 
= 25 and 6=1. 



approximate formula as 
P 



„,vro X[t) s 1 / g e-l°l g 2 (" +1 > cos (2|g| ^TTt) 



n=0 

-r „ oss t 



2(n + l) 



(11) 



Eq. is suggested as a quantum counterpart of the 
semiclassical solution, Eq. (JJJ, on the basis of compar- 
ison of the semiclassical and quantum solutions for the 
case 6 = 0, Eq. JHJ and Eq. JTUJ|. With T guess accu- 
rately guessed, Eq. ( fTTf l is expected to closely approxi- 
mate pte act (t): as l° n g as * is n °t extremely large. At very 
large times, Eq. (TTJ) predicts that p^P rox (t -► oo) -> ±, 
whereas the exact solution should yield p e e x e act (t — > oo) — > 
because of the irreversible "leak" into the lowest energy 
state 0) (see Fig- EJ - A similar type of leak was dis- 
cussed in an investigation of the collapse and revival in 
the absence of dissipation [l7j . 

The best guessed value for T guess on the basis of 
the quantum argument given above is T guess = ^T. 
On the other hand, the semiclassical theory predicts 
r guess = fr. The approximate solutions p q ^ ant {t) and 
plt ml (t), which are obtained by substituting T guess = 
and T guess = fr, respectively, into Eq. l[TTJl. are plot- 
ted in Fig. |3f a) along with the exact numerical solution 
Pte act (t)- We note that, on the scale with which this fig- 
ure is drawn, yofe° n *(i) is indistinguishable from p e e % aa '(t) , 
whereas p 8 e e e ml, {t) clearly exhibits a faster damping than 
pf e ac \t) and p q ™ n {t). Since the effect of different damp- 
ing rates shows up most clearly in the revival region, we 
show in Fig. &b) an expanded view of the three curves, 
P e ee act (t), P s et mi {t) and p%t ant {t), in the first and second 
revival regions. Even in this expanded figure, the two 
curves p e e xact (t) and pl^ ant (t) are indistinguishable. 

As a further comparison we plot in Fig. 01 the dif- 
ferences S semi (t) = p e e xact {t) - pf e ml {t) and 8 quant {t) = 
Pee ac *(*) - Ple ant (i)- while 5 semi (t) reaches a value as 
large as 3 x 10~ 2 , S quant (t) remains within 1CT 3 . We 
conclude therefore that the correct damping rate for the 
case b = 1 is in disagreement with the semiclassical 
prediction of fr. 

In conclusion we have shown that the Rabi oscilla- 
tion exhibited by a two-level atom interacting with a 



monochromatic field damps out with the rate of -j, re- 
gardless of the branching ratio b of the spontaneous decay 
into the ground state, where T represents the total decay 
rate of the excited state. This is in contradiction to the 
semiclassical prediction that the damping rate increases 

The reason 



with b from ^ at b = to 



|r at b = 1. 

for the constant damping rate lies in the openness of the 
quantum transition between the two states \g, n+1) and 
|e, n). The oversimplified structure of the transitions em- 
ployed in the semiclassical theory, as depicted in Fig. ^ 
cannot correctly describe this openness which is clearly 
visible in the full quantum mechanical structure shown 
in Fig. |21 The effect of a smaller quantum damping rate 
compared with the corresponding semiclassical damping 
rate shows up most clearly in the amplitude of the revived 
Rabi oscillations. 
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